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ABSTRACT 



By means of N-body simulations we investigate the impact of minor mergers on the angular momentum and dynamical properties of 
the merger remnant. Our simulations cover a range of initial orbital characteristics and gas-to- stellar mass fractions (from to 20%), 
and include star formation and supernova feedback. We confirm and extend previous results by showing that the specific angular 
momentum of the stellar component always decreases independently of the orbital parameters or morphology of the satellite, and that 
the decrease in the rotation velocity of the primary galaxy is accompanied by a change in the anisotropy of the orbits. However, the 
decrease aff'ects only the old stellar population, and not the new population formed from gas during the merging process. This means 
that the merging process induces an increasing diff'erence in the rotational support of the old and young stellar components, with the 
old one lagging with respect to the new. Even if our models are not intended specifically to reproduce the Milky Way and its accretion 
history, we find that, under certain conditions, the modeled rotational lag found is compatible with that observed in the Milky Way 
disk, thus indicating that minor mergers can be a viable way to produce it. The lag can increase with the vertical distance from the 
disk midplane, but only if the satellite is accreted along a direct orbit, and in all cases the main contribution to the lag comes from stars 
originally in the primary disk rather than from stars in the satellite galaxy. We also discuss the possibility of creating counter-rotating 
stars in the remnant disk, their fraction as a function of the vertical distance from the galaxy midplane, and the cumulative eff'ect of 
multiple mergers on their creation. 

Key words, galaxies: interaction - galaxies: formation - galaxies: evolution - galaxies: kinematics and dynamics 



1. Introduction 

Rotationally supported disks account for only a small fraction 
of the mass in the local Universe, but they contain most of 
the angular momentum (hereafter AM). The way disks acquire 
and redistribute their AM represents one of the most challeng- 
ing problems for models of galaxy formation and evolution. 
According to the current cosmological paradigm, baryons and 
dark halos in galaxies acquire their spin through tidal torques 
exerted by adjacent structures at early times. This AM is then 
redistributed among the different galaxy components through a 
number of internal and external processes as the galaxy evolves. 
Among the internal processes, bars, lopsidedness, spiral patterns 
and other coherent structures are efficient in redi stributing AM 
in ga l axies, as many studie s have shown (e.g 
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usually destroyed by the strong energy and AM redistri 
bution taking place during the interaction 
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). These stel- 



lar asymmetries can be stimulated or strengthened by external 
processes, such as accre tion of a few M^ yr~^ o f gas from cos- 
mological filaments (seelBournaud et al" 2()05allbb or tidal inter- 
actions and mergers (I Jog & Mavbhatel2006l:lMapelli et"^l2008l: 
iReichard et al .112009b . In particular, during an interaction orbital 
AM is converted into internal rotation, in an outside-in manner: 
the components which first interact are the most extended ones, 
while the more tightly bound components experience stro ng tidal 
eflTects only in the final phases of t he merging process (iBarnesI 
ll992l:lDi Matteo et al.ll2008bll2009l) . 

Many studies have shown that major mergers have a 
catastrophic impact on the ordered motion of the pre- 
existing galaxies. If the progenitors have disks, they are 
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Bendo & Barnes "2000; "Naab & Burkert* l2003l: |Bpurnaud elaD 
2005c; Jesseit et al . 2009), unle ss pecu liar orbital con figurations 
are chosen dPuerari & Pfennigej200ll:ICrocker et al.ll2009l) . The 
fraction of gas present in the progenitor disks can also influ- 
ence the morpholog y and kinematics of the final remnant (e.g., 
iHopkins et al ]|2009|), but this depends as well o n the gas physics 
implemented in the models (e.g.. iBournaud et al. 2011). In the 
case of pressure-supported progenitors, in turn, the tidal torques 
exerted by the companion can be strong enough to produce 
high rotational support (v/cr > 1) at large radii , even in merger 
remnants having an elliptical-like morphology (iDi Matteo et al.l 
120091) . More attention has been given to the study of the im- 
pact of AM redistribution in major mergers than in minor merg- 
ers (with mass ratios <0.1). This despite the fact that minor 
mergers are expected to be mu ch more common than major 
mergers (iFakhouri & Mai l2008l) and that many traces of on- 
goin g or past inter actions are visible both in the Milky Way 



(see iKlement l2010l for a recent review), our neighbor g alaxy 
Andromeda (llbata et al.l l200ll: iMcConna chie et al.' "2009) and 



other galaxies in the l ocal U niverse dMartinez-Delgado et al.l 
I2OIOI) . iDi Matteo et al.1 (l2010h pointed out that the AM redistri- 
bution during minor mergers also have an impact on the kinemat- 
ics of stellar disks and may explain the distribution of the orbital 
eccentricities of stars in the solar neighborhood. Understanding 
how AM is redistributed during such episodes is fundamental 
to the understanding of how disks can be maintained, how their 
kinematics can be aflfected, and what the signatures are of these 
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Table 1. Parameters for the initial models of the bulge, halo and 
stellar and gaseous disks. 





gso 


gSa 


gSb 


dEO 


dSO 


dSa 


dSb 


Mb [2.3 X 10^ Mo] 


10. 


10. 


5. 


7. 


1 


1 


0.5 


Mh [2.3 X 10^ Mo] 


50. 


50. 


75. 


3. 


5 


5 


7.5 


M, [2.3 X 10^ Mo] 


40. 


40. 


20. 





4 


4 


2 


Mgas/M, 




0.1 


0.2 






0.1 


0.2 


re [kpc] 


2. 


2. 


1. 


1.3 


0.6 


0.6 


0.3 


Th [kpc] 


10. 


10. 


12. 


2.2 


3.2 


3.2 


3.8 


a* [kpc] 


4. 


4. 


5. 




1.3 


1.3 


1.6 


[kpc] 


0.5 


0.5 


0.5 




0.16 


0.16 


0.16 


«gas [kpc] 




5. 


6. 






1.6 


1.9 


hgas [kpc] 




0.2 


0.2 






0.06 


0.06 



processes on the dynamical properties of different stellar popu- 
lations. 

This paper is the second of a series where we study, by means 
of numerical simulations, t he impact of min or mergers on AIM 
redistribution in galaxies. In lOu et aD (1201 Oh (hereafter Paper I), 
we investigated the impact of dissipationless minor mergers on 
disk galaxies, showing in particular that the initially non-rotating 
dark matter halo of the primary galaxy always gains AIM and 
that the specific AIM of the stellar component always decreases. 
We also showed that this decrease in AIM is accompanied by a 
change in stellar velocity anisotropy as the stellar orbits become 
less tangentially dominated as the merger advances. In this pa- 
per, we aim to advance this analysis by studying simulations of 
dissipative minor mergers, exploring a range of gas fractions and 
morphological parameters for the primary galaxy and the satel- 
lite. Star formation and feedback from supernovae explosions 
are included in the models, and we are able to trace the AIM re- 
distribution of all the galaxy components: dark matter, gas, old 
stars (i.e., those already in the galaxies before the interaction 
starts) and new stars (i.e., those formed from the gas during the 
interaction). In particular, we aim to understand if dissipative mi- 
nor mergers still slow down the stellar disk of the primary galaxy 
and if stellar populations of diff'erent ages show a diff'erent AIM 
content and diff'erent dynamical properties in the final (i.e., post- 
merger) disk. 

The paper is organized as follows: the numerical code, the 
initial galaxy models and orbital conditions adopted for the runs 
are described in §[2l Section [3] presents the main results, in par- 
ticular how the AIM content of gas and the old and new stellar 
populations is aff'ected by a single and by two consecutive (§ 13. lb 
minor mergers, and how the rotational lag induced by this redis- 
tribution aff'ects the old stellar component (§ 13.21) . The contribu- 
tion of accreted stars to the rotational lag, and the possibility to 
distinguish them on the basis of their AIM content, are discussed 
in § I3.3[ Section [T4l discusses the eff'ect of single and two con- 
secutive retrograde mergers on the fraction of counter-rotating 
stars. Finally, §|4]presents a discussion of the results found and 
the main conclusions. 



2. Models 

All 121 simulations described in this pap er are part of the 
GalMer projecQ and were fully described in IChilingarian et all 
(120101) . Here we recall the main characteristics of the adopted 
galaxy models and orbital parameters, as well as those of the 
numerical code employed to run the simulations. 



http://galmer. obspm.fr 



In this paper, we study the interaction and successive merger 
of a satellite galaxy with a ten times more massive primary disk 
galaxy. Both the primary and the satellite consist of a spheri- 

" cal non-rotating dark halo and a central bulge, both modeled by 

" Plummer spheres of masses IMh and 1Mb , respectively, and core 
radii rn and re, a stellar and an optional gas disk, represented by 
IMyamoto-Nagai density profiles of masses IM., and IMgas, disk 
scale lengths a* and agas, and scale heights h=, and hgas, respec- 
tively. Both the primary and the satellite galaxies span a range 
of morphologies, with a variety of bulge to disk ratios, and gas- 
to-stellar mass fractions (from fgas = to fgas = 20%). IMoving 
from early to late type systems, we refer to these models respec- 
tively with the nomenclature gSO, gSa and gSb for the primary 

" galaxy, and dSO, dSa and dSb for the satellite. For the satellites, 
we also consider a simple spheroid dominated model, consisting 
only of a dark halo and a stellar bulge, without any disk compo- 
nent - we refer to this model as dEO. All the parameters of the 
galaxies described above are given in Tabled and the number of 
particles adopted to describe the diff'erent galaxy components in 

Table[2l 

As we have done in lOu et all (l2010l) , all galaxy models were 
first evolved in isolation for 1 Gyr before the interaction starts, to 
let the initial system reach a stable configuration. Once relaxed, 
the two galaxies were placed at an initial distance of 100 kpc, 
wit h a variety of re lative velocities, to simulate diff'erent orbits 
(see IChilingarian e t al. 2010, Table 9, for the orbital initial con- 
ditions). To study the eff'ect of multiple mergers, we also ran 
some simulations in which the primary galaxy accretes consec- 
utively two or three satellites over a period of 3-5 Gyr. We also 
compared the dynamical properties of disks after minor mergers 
to those heated by internal processes, such as sc attering by mas- 
sive c lumps formed in an initially unstable disk (iBournaud et al.l 

l2009h . For this, we analyzed a si mulation of a gas-rich, unstable 
disk galaxy from lDi IMatteo et al . ( 2008a). It uses a total number 
of 120,000 particles, equally distributed among gas, stars and 
dark matter, and has the same parameters as the gSb model (see 
Table[T]) except for the gas mass fraction, which is initially much 
higher, 50%. 

All th e simulations were mn usi ng the Tree-SPH code de- 
scribed in lSemelin & CombesI (l2002h . Gas is treated as isother- 
mal, at a fixed temperature of T=10^ K. Prescriptions for star 
formation and feedback from supernovae explosions are also in- 
cluded. The rate of star formation is governed locally by a vol- 
ume density Schmidt law IMgas ^ P^Mgas with a=l.5, which is 
implemented in the code by means of a hybrid particle ap proach 
(^see llMihos & HernquistI [l994l: IChilingarian et al.ll2010( for de- 
tails). At each time step, a hybrid particle is characterized by two 
mass values: the first is its total mass M/ which stays constant 
during the simulation and is used to calculate the gravitational 
force, and the second is given by the gas content of the particle, 
Mi,gas, which changes with time according to the Schmidt law 
and is used to evaluate hydrodynamical quantities. At each time 
step the mass of the new stars, i.e. those formed since the be- 
ginning of the simulation, is given by IMi -IMgas. If Mi^gas drops 
below 5% of the initial value, the "hybrid" particle is considered 
to be totally converted into a star-like particle and its small resid- 
ual amount of gas spread out over its (hybrid) neighbors. Note 
that in the hybrid particle scheme, new stars follow the gas kine- 
matics until the moment when the hybrid particle they reside in 
is completely converted into stars. This can represent a limitation 
of our approach, in the sense that it does not allow us to prop- 
erly follow the heating with time of the newly formed stars by 
secular processes. Finally, for the evaluation of the gravitational 
forces a softening length 6=200 pc is employed. The equations 
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Table 2. Particle numbers for primary disk galaxies and satel- 
lites. 
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dEO 


dSO 


dSa 
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N 

Nstar 

Ndm 


32,000 
16,000 


32,000 
16,000 


8,000 
24,000 
16,000 


16,000 
16,000 
16,000 



of motion are integrated using a leap-frog algorithm, with a fixed 
time step of 0.5 Myr. With these choices, the relative error in the 
conservation of the total energy is close to 10"^ per time step. 

The nomenclature adopted for our simulations consists of a 
string of 13 characters: the first three (gSO, gSa or gSb) for the 
primary galaxy and the following three for the satellite galaxy 
(dEO, dSO, dS a or dSb), followed by t he encounter identifica- 
tion string (see lChilingarian et"ani2010l Table 9) and the orien- 
tation of the disk of the primary galaxy with respect to the orbital 
plane (33 or 60 degrees). Multiple mergers are indicated by two 
more characters - tn, where n is the number of satellites accreted 
by the primary galaxy. Eor example, the string gSadSat202dir33 
refers to a gSa galaxy accreting two dwarf dSa galaxies, whose 
initial orbital para meters are those correspo nding to the string 
02dir in Table 9 of IChilingarian et al.l (l2010h . 

3. Results 

In Paper I, we showed that during dissipationless minor mergers 
the orbital AM is redistributed into internal AM of the interact- 
ing galaxies in an outside-in manner: the initially non-rotating 
dark halos of primary and satellite start to acquire AM first, 
just after the first pericenter passage between the two galaxies, 
while the strongest changes in the primary stellar disk take place 
mostly in the final phase of the collision, when the satellite is 
close to merging. All our simulations also showed that the stel- 
lar disk loses part of its initial AM. The aim of this section is to 
investigate if these results are also valid for dissipative mergers, 
if diff'erences can be found in the rotational support of the dif- 
ferent stellar populations in the primary disk once the merger is 
completed and what their main kinematic characteristics are. 

3.1. Gas, old stars, new stars and the evolution of their 
angular momenta 

As discussed in §[2l our model galaxies contain both "old stars", 
which were already present before the start of the interaction, 
and "new stars" formed from gas during the interaction. We will 
follow the nomenclature of old and new stars throughout the text. 
An example of the spatial distribution and morphological evolu- 
tion of these two stellar components, as well as the dissipative 
gas component, is shown in Eig[T] Before the first close passage 
of the satellite, which occurs at 0.5 Gyr, gas and new stars are 
distributed in a morphologically thin and dynamically cold disk, 
and old stars in a thicker and hotter component and in a central 
bulge. At 0.65 Gyr, the primary disk shows signs of a tidal per- 
turbation: some gas and new stars are clearly displaced from the 
midplane of the primary, and the old stellar disk starts to thicken 
even further. Note also that the satellite morphology has already 
been strongly modified at this time. Signs of tidal perturbation of 



the primary galaxy are visible all along the merging sequence af- 
ter the first close passage of the satellite. At 2.6 Gyr, or more than 
1 Gyr after the merger ended, its morphology has clearly been 
profoundly modified: a thic k disk made of old stars has been 
formed (see lOu et al.|[201ll for a description of its morphologi- 
cal properties) whereas gas and new stars are still distributed in a 
thin disk component with some perturbations at large radii, such 
as a warp and a tidal tail. 

During the interaction, the eff'ects of dynamical friction and 
tidal forces redistribute the AM in the system: both the primary 
and the satellite galaxy acquire part of the orbital AM, and, as 
already found in Paper I, the most extended regions experience 
these changes first. In particular, the impact of a minor merger on 
the AM content of the primary galaxy is shown in Eig.[2l where 
the evolution with time of the specific AM of old stars, gas, new 
stars and dark matter is shown for diff'erent regions of the galaxy. 
The figure shows clearly that: 

- the outermost regions of the initially non-rotating dark halo 
are the first to acquire part of the orbital AM, already after 
the first pericenter passage, at ^ = 0.5 Gyr, while inside 3Rd 
(the initial disk scale length 7?^ is 4.8, 4.8 and 5.2 kpc for SO, 
Sa and Sb type galaxies, respectively) tend to increase their 
AM only during the final phases of the merger; 

- the AM of the baryonic components is also aff'ected by the 
interaction. The old stars experience a slowing down of their 
rotation at all radii, especially outside O.SRd. A less pro- 
nounced decrease of the AM is visible in the gas (and new 
stars) between and 3Rd. For example, in the merger case 
shown in Eig. [2l the fractional decrease in the specific AM 
of old stars is A/// = 0.19 at 2-3Rd whereas A/// = 0.09 for 
gas (and new stars) in the same region. While the outermost 
regions gain AM in the innermost ones the specific AM re- 
mains unchanged. 

The slowing down with time of the AM content of old stars 
at a given radius is not only due to disk heating, but also an 
eff'ect of the merger process (as has been previously shown by 
Quillen et al. 2009). To demonstrate this eff'ect, in Eig. [3] we 
show the total AM of old stars, gas and new stars (for dissi- 
pative mergers), and dark matter of the primary galaxy for three 
simulations with diff'erent gas fractions in the primary disk (0, 
0.1 and 0.2). In all cases, one can see a substantial decrease of 
the total AM of old stars (AL/L = 0.13 for a gSbdSb direct 
merger, 0.2 for a gSadSa merger and 0.36 for a dissipationless 
merger). Also the gas (and new stars) tend to lose a fraction of 
their AM. At the same time, the dark halo of the primary galaxy 
acquires part of the AM, which is already apparent after the first 
passage of the satellite, around t = 0.5 Gyr in the simulations. 
These substantial changes in AM content of the diff'erent com- 
ponents are clearly an eff'ect of the minor merger itself. This can 
be seen when comparing the evolution of the merging galaxies 
with that of the corresponding galaxies which were evolved in 
isolation (right panels of Eig. O. Secular processes, which dom- 
inate the evolution of the isolated galaxy simulations, cause less 
pronounced changes in the AM, with the decrease in AL/L of the 
old stellar disk being only about 20-25% of what is estimated for 
the minor merger simulations. This demonstrates that disk heat- 
ing due to secular processes, which increases the radial stellar 
velocity dispersions, cannot be the only mechanism responsible 
for the slowing down of the disks observed in the minor merger 
simulations. 

Eigs.[2]and [3]suggest that it should be possible to find a dif- 
ference in the AM content of old and new stars in the final, post- 
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Fig. 1. From left to right: evolution with time of gas (upper panels) and stars (middle and lower panels) during a minor merger 
between a gSb and a dSb galaxy on a direct orbit (gSbdSb01dir33). The stellar component is separated into new stars (formed 
during the interaction, middle row of panels), and old stars (those already present in the galaxies before the interaction, lower 
panels). In this projection the disk of the primary galaxy is viewed edge-on. 
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Fig. 2. Evolution with time of the specific AM, /, of a gSb pri- 
mary galaxy during the minor merger whose evolution is shown 
in FigdJ The specific AM of old stars (left panel), gas and new 
stars (middle panel) and dark matter (right panel) is shown, for 
five diff'erent radial zones in the galaxy, see the legend. Rd indi- 
cates the disk scale length. The specific AM is in units of 100 
kpc km s"^. 



merger disk, as is indeed the case, see Fig.|4l comparing the spe- 
cific AM as a function of radius for gas and new stars (left panel) 
and old stars (middle panel), one can see that while the AM of 
the new stellar component is unchanged after the merger, that of 
the old stars decreases at all radii (A/// = 0.23 at Rd). This is due 
to the fact that, during the merging process, old stars are heated 
as shown by the increasing radial and vertical velocity disper- 
sion (Fig. [5]), which then leads to a slowing down of the stellar 
disk, as we have argued in paper I. However, gas can dissipate its 
energy, and thus can preserve mainly tangential motion and keep 
its AM unaff'ected, whereas stars newly formed from this gas re- 
tain the AM that was acquired by the gas. This means that at a 
given radius the old and the new stellar populations have a diff'er- 
ent specific AM and thus diff'erent amounts of rotational support. 
As discussed in Paper I for dissipationless mergers, the decrease 
in the rotation speed of the old stellar disk of the primary galaxy 
is accompanied by a change in the distribution of the types of 
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Fig. 3. Left panels: Evolution with time of the total AM of gas + 
new stars (dashed line), old stars (dotted line) and the dark matter 
halo (dashed-dotted line) in the region at r<20 kpc from the cen- 
ter of the primary galaxy during a minor merger with a satellite 
galaxy on a direct orbit. From top to bottom, the gas mass frac- 
tion of both primary galaxy and merging satellite are 20%, 10% 
and 0, respectively. The AM is in units of 2.2xlO^^Mo kpc km 
s"^ Right panels: The evolution of the AM of the corresponding 
components (gas -h new stars, old stars and dark matter halo) of 
isolated galaxies with the same initial structure as the primary 
galaxies in the minor merger simulations. 



Qu et aL: Minor mergers and their impact on the kinematics of old and young stellar populations in disk galaxies 



5 



gas 


old stars 






old stars 




_ new stars 















before merger/^ 




-2 












-4 


\ post- 


merger _ 








-6 


V„„^^before 


merger _ 




/^!^^ost - m e rg e r _ 




-8 






1 2 3 


1 2 3 






1 2 


3 










r/R. 





t = 0.05Gyr 



gSbdSb01dir33 
0.5Gyr 



1000 
100 
10 




before merger 



2 



2 



Fig. 5. Left panel: Ratio of tangential and radial velocity dis- 
persion as a function of scaled radius r/Rd for old stars before 
(black line) and after (red line) the minor merger shown in Fig[T] 
Right panel: Same as the left panel but for the ratio of vertical 
and radial velocity dispersion. 



Fig. 4. From left to right: Specific AM as a function of radius 
normalized by the disk scale length, r/7?j, of the primary galaxy, 
for gas and new stars {left panel) and old stars {middle panel) 
during the minor merger simulation whose evolution is shown 
in Figd] In each panel the black line shows the initial AM and 
the red line shows the final AM. The specific AM is in units of 
100 kpc km s"^ Right panel: The anisotropy parameter as a 
function of r/Rd for old stars in the same merger. Both the initial 
(black line) and final (red line) p are shown. 
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Fig. 6. Upper panels: Relative change of the specific AM with 
respect to the initial value. Mint /lint = {LtO) - LtO = ^))llint{t = 
0), of old stars of a gSb galaxy as a function of initial radius 
rini (scaled by Rj) during a minor merger with a dSb galaxy on 
a direct orbit (gSbdSb01dir33). Also shown is the average of 
Mint I lint as a function of n^^i/Rd in each epoch (overlaid solid 
line). The vertical bar on the left shows the color scale indicating 
the number of old stars at certain MintlUnt values. Lower panels: 
The distribution of Mintlhnt of old stars of an isolated galaxy 
with the same initial structure as the primary galaxy in the minor 
merger case. 



stellar orbits: the radial component of the velocity dispersion be- 
comes increasingly important during the merger, thus increasing 

2 

the anisotropy parameter, /3 = 1- ^rr . from its initially negative 
value. This is also the case for dissipative minor mergers (Fig.©. 
Independent of the amount of gas present in the primary disk or 
in the satellite galaxy, minor mergers always result in a slowing 
down of the old stellar disk. The fractional decrease in the spe- 
cific AM at Rd is A/// = 0.23 for a gSb model. A/// = 0.17 for a 
gSa and A/// = 0.29 for a gSO. 

Note that the slowing down of the old stellar disk during 
the merger is also visible in the time evolution of the AM of 
those stars that were initially (at t=0) present in a given radial 
bin, rather than in that of stars currently present at a given ra- 
dius (Fig.O. At any given initial radius, r/„/, the variation of the 
AM is much stronger when the stellar disk undergoes a minor 
merger than when it evolves in isolation. In the minor merger 
simulations, the average value of A/// decreases with time, espe- 
cially for stars whose initial radii exceed Rd. The largest changes 
take place in the outer disk where the eff'ect of the perturbations 
due to the satellite i s the greatest (as already has been shown by 
iQuillen et al.ll2009l:lBird et al.ll201lh . 

3.1 .1 . The impact of multiple mergers 

If a single minor merger causes the slowing down of the old 
stellar component of a galaxy disk, what would be the eff'ect of 
a subsequent minor merger on it? To answer this question, we 



ran ten simulations in which the primary galaxy consecutively 
accretes two ten times less massive satellites over a period of 3- 
5 Gyr. The second minor merger produces a further decrease in 
the specific AM of the old stellar population (Fig.l?]). After the 
first merger, the specific AM decreases at all radii, accompanied 
by an increase of the ft parameter. In the second merger event, 
the AM content of the old disk stellar component is further de- 
creased and this is accompanied by an increase in the velocity 
anisotropy parameter JS at all radii in the disk. Furthermore, we 
ran one simulation in which the primary galaxy consecutively 
accretes three ten times less massive satellites over a period of 
5 Gyr. After having accreted two satellites the stellar disk of the 
primary galaxy still reacted dynamically to the third satellite ac- 
cretion - its specific AM decreased further, its beta parameter 
increased and the slowing down of the old stellar disk did not 
saturate but increased (see Fig. (T]). While old stars slow down 
with time, whereas stars formed during the repeated mergers 
tend to show little change in their AM (Fig. (8]), the diff'erence 
in AM content between stars born at diff'erent times increases 
during repeated mergers with low mass satellite galaxies. This 
is only partially due to the hybrid method used to implement 
star formation in our simulations, in which the new stars are 
continuously formed from dissipative gas which keeps the AM 
mostly unchanged. Even though secular disk heating processes 
also have an eff'ect on heating newly formed stars, this eff'ect is 
significantly lower than that induced by minor mergers, as al- 
ready shown in Fig. [3] 
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Fig. 7. Left panel: Specific AM as a function of scaled radius 
rjRd for the old stars of a gSa galaxy undergoing three consecu- 
tive minor mergers. The specific AM is shown after the first (dot- 
ted line), the second (dashed-dotted line), and the third (dashed 
line) merger, and the initial AM is shown for comparison (solid 
line). The specific AM is in units of 100 kpc km s"^ Right panel: 
Anisotropy parameter, yS, as a function of scaled radius, for the 
old stars of the gSa galaxy whose specific AM is shown in the 
left panel. 
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Fig. 8. Same as Fig. [7] but for newly formed stars in the gSa 
galaxy. 

3.2. Rotational lag 

As the slowing down is only aff'ecting the old stellar population, 
we expect to see an increasing diff'erence between the AM con- 
tent of the old and new stellar components as mergers proceed. 
This is indeed the case, as shown in Figs. l9l and [TOl where the 
distribution of the specific AM, /, and radial and tangential ve- 
locities ivr and Vt, respectively) of stars between radii of 2 and 
?>Rd are shown after a single minor merger and after two consec- 
utive mergers. We chose this radial region as it should be rep- 
resentative of the solar neighborhood. To compare the eff'ects of 
consecutive mergers to those of secular evolution, we show in 
Fig. [TT] the same distributions for three galaxy simulations that 
do not undergo a merger but are evolved in isolation: for gSa and 
gSb model galaxies, and for a gas-rich, unstable clumpy disk. In 
this plot we show the newly formed stars as well as the entire 
old star population and the contribution of old stars at heights 
I z |<1 kpc, in order to distinguish old stars distributed in a thin 
disk from the total. We emphasize that this is a very simple way 
to distinguish between a thin, dynamically cold stellar compo- 
nent and a thick and hotter stellar component, which is based on 
the analysis of the vertical pro perties of thick disks formed in 
minor mergers (lOu et al.ll201ll) . 

The results of this study show that: 

- After a single merger the thin stellar disk (i.e., at | z |<1 kpc) 
consists of two diff'erent stellar components, a young and an 
old population, the latter with a lower angular momentum 
than the former. In the region between 2 and 37?j, for a gS- 
bdSb direct merger, where both galaxies initially have a gas 
fraction of 20%, the old stars in the thin disk have an average 
specific AM which is about 10% lower than the new stars; 



- If the whole old stellar population is considered, its aver- 
age specific AM is about 7% lower than that of old stars at 
I z |<1 kpc, indicating that the stars in the thick old disk 
(I z |>1 kpc) contribute to decrease the specific AM and are 
thus less rotationally supported than the old stars in the thin 
disk; 

- This diff'erence in specific AM content of the stellar popu- 
lations is reflected in difl'erent tangential velocities of their 
stars in the disk, with the tangential support of stellar orbits 
decreasing from the thin new disk, through the thin old disk, 
to the thick old disk; 

- For the same gas fraction in the progenitor disks there is 
a difl'erence between direct and retrograde mergers, in the 
sense that the specific AM content of old stars is lower if the 
satellite orbit is retrograde. This lower AM content can be 
explained by the presence of counter-rotating material (neg- 
ative / ) which is found in all retrograde encounters analyzed 
(see § l3.4l for a detailed discussion); 

- If one defines the rotational lag of old stars as viag = - 

new star whcrc V* is the tangential velocity of stars in the 
thin or thick old disk, and v^w stars is that of new stars, one 
can see that it depends on the amount of gas present in the 
progenitor disks: for old thin and old thick disk stars, the 
higher the gas fraction in the progenitor disk, the lower the 
rotational lag in the remnant. 

We refer the reader to Table [3] for a complete summary of the 
results found in Figs. [9l[TQl and [TTl 

The diff'erence between the specific AM content of old and 
new stars becomes still more pronounced after a second, consec- 
utive minor merger, leading to a further increase in the rotational 
lag of both the old thin and thick disk stars, as shown in Fig.fTOl 
Given that minor mergers have a cumulative impact on heating 
and slowing down stellar disks (see also Fig. [7]), two mergers of 
mass ratio 1:10 are expected to result in a diff'erence between 
the specific AM content of old and new stars that is comparable 
with that of one merger of mass ratio 1:5. In other words, more 
massive satellites can have a larger impact on slowing down the 
stellar disk and produce a larger rotational lag betwee n the old 
and n ew stellar populations than less massive satellites (lOu et al.l 
I2OIOI) . The substantial diff'erence between the specific AM con- 
tent of the thin and thick disk stellar populations cannot be due 
to secular evolution processes alone, which produce a consid- 
erably smaller variation in both / and Vt in 3 Gyr of evolution 
(Fig. [TT]). Note also that secular evolution produces a narrower 
distribution of /, Vr and v^, thus suggesting that also mixing and 
radial migration due to secular processes may be less eff'ective 
than in minor mergers. Whereas scattering of stars by massive 
clumps can induce a rotational lag (in old thin and thick disk 
stars), the overall impact of scattering by mass concentrations is 
(moderately) lower than that produced by a direct 1:10 merger 
on a galaxy with an initial gas fraction of 20%. 

3.3. Accreted and heated disk stars 
3.3.1 . Contribution to the rotational lag 

In the previous section we saw that old stars in the thin and thick 
disk of a minor merger remnant lag with respect to the new stars, 
and that this lag is higher for stars further from the galaxy mid- 
plane (I z |>1 kpc). But what is this rotational lag due to - is 
it mostly associated to stars originally in the satellite or in the 
primary disk, and how does it depend on the orbital parameters? 

To answer these questions, we analyzed the tangential veloc- 
ity Vt of old stars at radii between 2 and 3Rd in four diff'erent 



Qu et aL: Minor mergers and their impact on the kinematics of old and young stellar populations in disk galaxies 



7 



gS0dS001dir33 

old stars 
IzKlkpc 

old stars all 



-40 



-20 




' 1 ' 

gSadSa01dir33 

new stars 

old stars 
IzKlkpc 

old stars all 




-40 



-20 



20 



40 



gSbdSb01dir33 

new stars 

old stars 
IzKlkpc 

old stars all 




-40 



-20 



20 



40 



gS0dS001ret33 

old stars 
IzKlkpc 

old stars all 



■/////A 



-40 



-20 




gSadSa01ret33 

new stars 

old stars 
IzKlkpc 

old stars all 




-40 



-20 



20 



40 



' 1 ' 1 
L gSbdSb01ret33 


' II 1 ' - 

II 1 : 
II 1 


I new stars 




r old stars 


II 1 J 


: IzKlkpc 


r old stars all 




1 





-40 



-20 



20 



40 













Vt 



Fig. 9. Histograms of the specific AM, / (in units of 100 kpc km s"^), and the radial (v^) and tangential (v^) velocities (in units of 
100 km s"^) of stars between radii of 2-37?j in a number of minor merger remnants. The top three rows show direct minor mergers, 
with disk gas mass fractions increasing from to 0.2. Indicated in each panel are the entire old stellar population (gray), new stars 
(blue) and old stars with | z |<1 kpc (red). The vertical dashed lines show the average values of /, and for each of these three 
components. The three lowest rows are similar, but for retrograde mergers. 
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Fig. 10. Same as Fig.|9l but for two consecutive minor mergers. 
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Fig. 11. Same as Fig. [51 but for galaxies that were evolved in isolation for 3 Gyr: a gSa (upper panels) with a gas mass fraction 
/gas = 0.1 and a gSb (middle panels) with /gas = 0.2, as well as a galaxy with a much higher gas fraction of 0.5 which went through 
an unstable clumpy phase (bottom panels). In the latter case the parameters are also shown after 3 Gyr, when the clumps have 
dissolved through tidal effects or spiraled into the galaxy center through dynamical friction and interaction. 



regions above the disk midplane: | z |<1 kpc, 1<| z |<3 kpc, 
3<| z |<5 kpc and 5<| z |<10 kpc. Results from some represen- 
tative mergers are shown in Fig. [121 where for each encounter 
we distinguished the primary stars from those originally in the 
satellite. This figure shows some interesting trends. First of all. 



there is a clear diff'erence between direct and retrograde mergers. 
While in both cases the average Vt of old stars is lower than that 
of new stars (Table [3]), the behavior of Vt as a function of height 
z is significantly diff'erent. Direct mergers produce disks with a 
tangential velocity Vt which decreases with height, meaning that 
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Fig. 12. Distribution of tangential velocities Vt (in units of 100 km s"^) of old stars for six minor mergers. For each encounter a 
pair of panels is shown, with stars originally in the primary disk (left) and those accreted from the satellite (right). In each panel 
the histograms correspond to different heights above and below the galaxy midplane: | z |<1 kpc (black), 1<| z |<3 kpc (blue), 
3<| z |<5 kpc (green) and 5<| z |<10 kpc (orange). The average values of Vt in the different regions is also indicated. 
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Table 3. Kinematics of stars in the post-merger thin and thick old thick disks and of the new stars in disk galaxies heated by different 
physical processes. For old stars in the thin and the thick disks the rotational lag is defined with respect to the tangential velocity Vt 
of the new stars. 
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Fig. 13. Same as Fig.[T2]but for an isolated, very gas-rich galaxy 
(/gas = 0-5), after an unstable, clumpy phase. Note that in this 
case the fourth and highest region is missing, due to a lack of 
stars above/below 5 kpc from the galaxy midplane. 



the rotational lag increases with z, whereas in retrograde merg- 
ers Vt is constant with height up to z~5 kpc, and decreases only 
at greater heights. The satellite stars, which in each region con- 
stitute only a small percentage of the total stellar content, do not 
contribute significantly to the lag. IMoreover, if one compares the 
tangential velocities of satellite and primary stars for direct or- 
bits, one can find a variety of behaviors in the regions analyzed: 
in some cases the tangential velocities of satellite stars are al- 
ways smaller than those of stars from the primary (e.g. for orbit 
"gS0dS001dir33"), while in other cases the values are compara- 
ble, and satellite stars show even higher velocities than primary 
stars in the outer regions (e.g. orbit "gSadSa01dir33"). In all the 
cases, however, the strongest variations in Vt as a function of z 
are associated with stars originally in the primary rather than 
from the satellite. For retrograde orbits, the tangential velocities 
of satellite stars show only small variations with increasing z, 
as is the case for primary stars. We note also that a stellar thick 
disk formed in an unstable clumpy galaxy is characterized by 
tangential velocities whose variation with z is very similar to 
that produced in a direct encounter with /gas = 0.2 (compare, 
e.g. Fig. [13] with orbit "gSbdSb01dir33" in Fig.[l2]). 



formed in the primary thin disk and then heated by a minor 
merger from those originally in the satellite and then accreted. 
Their simulations showed (see their Fig. 6) that satellite stars in 
the height zone | z |<1 kpc have a typical AIM content which 
is constant with radius, while stars from the primary disk tend 
to have an AIM / oc r. Their simulations are gas-free and in- 
volve satellites more massive than the ones studied in this pa- 
per. Because of these diff'erences we decided to reinvestigate 
the validity of their method, exploring if the result is similar for 
larger mass ratios like ours (1:10) and for other satellite and pri- 
mary galaxy morphologies. Our results for the specific AIM of 
stars originally in the primary disk and in the satellite galaxy are 
shown in Figs.[T4]and [151 We divided the stars into two regions, 
one at I z |<1 kpc and another at | z |>1 kpc. The general result 
confirms that the diff'erence between the specific AIM of satellite 
stars and primary disk stars is significant in the outer disk re- 
gions: typically the diff'erence is 20% at radii of 1.5-27?j. But at 
these radii the fraction of satellite stars with U |<1 kpc is only a 
few percent of the total number of stars. This means that it can be 
extremely difficult to detect these stars in an observational sam- 
ple. We find, however, that the fraction of satellite stars is five to 
ten times larger at greater heights: at | z |>1 kpc, outside a radius 
of \.5-2Rd the specific AIM content of accreted stars can be sig- 
nificantly different from that of primary stars, and their fraction 
sufficiently higher (around 10%) to make the observational de- 
tection of this signature more likely. Note however that generally 
the satellite stars with the lowest AIM content are found very far 
from the galaxy midplane (| z |>5 kpc). We want to point out that 
there are also cases in which there is no discernible difference in 
the specific AIM of satellite and primary disk stars at any radius, 
even at the largest radii plotted in Figs. [141 and [15]). This is in 
agreement with what was found in Fig. [121 where primary and 
satellite stars show remarkably similar tangential velocities. 

Obviously, the specific AIM becomes a better discriminant 
if satellite orbits counter-rotate with respect to the primary disk 
- a clear sign of an external origin of stars. We find counter- 
rotation in all merger remnants which formed from satellites on 
retrograde orbits (some examples are shown in Fig.[T5l). We will 
discuss the main properties of counter-rotating stars further in 
the next section. 



3.3.2. How to distinguish satellite stars 

IVillalobos & Helmil (l2009l) recently proposed an interesting 
method to use the vertical AIM to discriminate between stars 



3.4. Counter-rotation 

Coun ter-rotating stellar disks are observed in external galaxies 
(see lYoachim & Dalcantonll2005l) and counter-rotation is also 
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Fig. 14. Specific AM, /, as a function of scaled radius for old stars in three direct merger remnants. For each encounter a pair of plots 
is shown, where stars are selected on their heights above (or below) the galaxy midplane: U |<1 kpc (upper panels) and | z |>1 kpc 
(lower panels). In each panel stars originally in the primary galaxy are shown as black points, stars originally in the satellite as blue 
points, and satellite stars with | z |>5 kpc as red points. The specific AM is in units of 100 kpc km s"^ . 
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Fig. 15. Same as Fig. [141 but for three retrograde mergers. 



found in the Milky Way halo (see ICaroUo et a02007l) . The pres- 
ence of stars (in the thick disk or in the halo) rotating in a direc- 
tion which is reverse of the galaxy main spin cannot be explained 
with secular evolution processes, and is a strong evidence for an 
external origin. 

It is thus of great interest to investigate the characteristics 
of counter-rotating stars in our merger remnants - specifically, 
their fraction and vertical distribution with radius. Examples are 
shown in Figs. [161 and [T7] for a number of dissipationless and 



dissipative retrograde mergers. In these plots four diff'erent re- 
gions are selected: (1) at | z |<1 kpc, (2) at 1<| z |<3 kpc, 
(3) at 3<| z |<5 kpc and (4) at 5<| z |<10 kpc. It is clear that 
counter-rotating stars are found at all values of z, from the galaxy 
midplane to 10 kpc above it, and that the probability to find 
counter-rotating stars increases with height: the fraction Net IN tot 
of counter-rotating stars to the total number of stars in the re- 
gion is <5% up to I z |<3 kpc, typically around 10% in the 
3 kpc<| z |<5 kpc and it can reach 30% in the region highest 
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Fig. 16. Fraction of counter-rotating stars (shown in blue) for four dissipationless retrograde mergers. Stars are divided into four 
regions, according to their height above (or below) the galaxy midplane: (1) at | z |<1 kpc, (2) at 1 kpc<| z |<3 kpc, (3) at 
3 kpc<| z |<5 kpc and (4) at 5 kpc<| z |<10 kpc. For each region the fraction of satellite stars to the total number of stars is also 
shown (gray). 
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Fig. 17. Same as FigO but for four dissipative retrograde mergers. 
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Fig. 18. Same as Fig[T6l but for a gSO galaxy after a single ret- 
rograde merger (left panel) and after two consecutive retrograde 
mergers (right panel) with a dEO satellite. 



above the disk. This trend is also due to the fact that the fraction 
of satellite stars A^^ to the total number of stars increases with z. 
In other words, it is the overall increase in the fraction of satellite 
stars with z which determines this trend, not simply a decrease 
in the number of disk stars as a function of height above the disk. 
Furthermore the fraction of counter-rotating stars depends also 
on the initial orbital inclination of the satellite: for a satellite on 
an orbit inclined by 60 degre es with respect to the prirn ary disk, 
this fraction is lower at all z (IVillalobos & Helmill2008h . 

If counter-rotating stars are a natural outcome of low orbital 
inclination in single minor mergers with satellites on retrograde 
orbits, would a second retrograde merger then increase the frac- 
tion of counter-rotating stars even further at all heights? In our 
models (Fig. [18]) the effect of a second retrograde merger with 
the same mass ratio is to increase the fraction of counter-rotating 
stars by a factor of about two, at all heights. We can conclude 
from this that repeated retrograde minor mergers have a cumu- 
lative effect on the number of counter-rotating stars in the thick 
disk and inner halo. 

3.5. Kinematics of the remnant disk 

The AM evolution of the baryonic component during merging 
processes discussed in the § 13. H is naturally reflected in an evo- 
lution of the line-of-sight velocity vios and the vioslcTy ratio. 
Minor mergers always produce a decrease in vios and an in- 



crease in the velocity dispersion cFy of the stellar component 
at all radii (Fig. [191). The remnant stellar disks usually have a 
smaller viosIcTv ratio than the initial values before interaction. 
Consecutive satellite accretions have cumulative effects on de- 
creasing the Vios and increasing the CTy, and thus deceasing the 
^los/o'v ratio. Nevertheless, smaller variations in both vios and 
(Ty can be found in minor mergers with gas, compared to gas- 
free merger cases, due to the fact that dissipative gas can help 
preserve disk rotation by forming new rotating disk stars during 
mergers. The decreases in vios and the vios/o-y ratio indicate that 
stellar disks become hotter and rotate more slowly after minor 
merger events, which is consistent with the AM loss and the in- 
crease of the anisotropy parameter (Figs. |4] and [7]). 

4. Discussion and Conclusions 

By means of N-body/SPH simulations, we studied the kinemat- 
ics of stars in galaxy remnants of 1 : 10 mass ratio minor mergers. 
The simulated interactions span a range in gas mass fractions, 
from to 0.2 in the primary and satellite galaxies, and a range of 

orbital parameters^ 

As shown in lOu et aP (l2010l) . minor mergers result in a re- 
distribution of orbital into internal angular momentum, which 
affects all galaxy components. In particular, old stars, i.e. those 
already in place before the interaction, always lose angular mo- 
mentum during the merging process. The decrease of the specific 
AM of old stars is accompanied by a redistribution of stellar 
orbits, as traced by the anisotropy parameter JS, which become 
increasingly radial. In minor mergers with gas in the disk of the 
primary galaxy we find a similar trend of old stars losing angular 
momentum and as a result their orbits becoming more radially 
dominated. However, when a new stellar component forms from 
gas present in the primary disk during the merger, its AM content 
is significantly different: the orbits tend to be more tangentially 
dominated, thus providing a higher rotational support. This dif- 
ferent behavior results in a final stellar disk with two different 
stellar populations with significantly different AM content. In 
particular, old stars always show a rotational lag with respect to 
the young stellar component. If one separates all stars into thin 
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disk stars (at heights | z |<1 kpc from the galaxy midplane) and 
thick disk stars (at | z |>1 kpc), three different components can 
be found, with different dynamical properties: (1) young stars 
in the thin disk, which are rotationally supported and show the 
highest values of v^, (2) old thin disk stars lagging with respect 
to the new stars and (3) old thick disk stars lagging with respect 
to both thin disk components. For a minor 1:10 merger, with a 
satellite accreted on a direct orbit and with an initial primary 
disk gas fraction of 0.2, the old stars in the thin disk have a ro- 
tational lag of about 20 km s"^ while the old stars in the thick 
disk have a velocity about 50 km s"^ lower than the young stellar 
component, both lag val ues being compa tible with the estimates 
for the Milky Way (see iGilmore et al.ll2002) . Multiple mergers 
can further reduce the tangential velocity of the old stellar com- 
ponents, while leaving that of the new stars mostly unchanged, 
thus resulting in a further increase in rotational lag with every 
successive accretion episode. As the two populations, old stars 
and new stars, have different tangential velocities, in a plot of Vt 
as function of age we expect to find a discontinuity at the time 
when the merger occurs. Of course, the newly formed stars will 
evolve in time, and if no other merger takes place they will be 
slowly heated by secular effects. However, as we have shown in 
Fig. [3 secular processes are much less effective in altering stellar 
kinematics than minor mergers. Therefore, in our opinion, a dis- 
continuity in the age-v^ plane (or the age-yS parameter) between 
the old and new stellar populations should still be visible, even if 
secular processes, and asymmetric drift in particular, contribute 
to the slow heating of the new stellar populations. Unfortunately, 
the hybrid particle method we adopted to implement star forma- 
tion in the simulations does not allow us to follow properly the 
heating of the newly formed stars with time. We do think, how- 
ever, that this discontinuity in the age-v^ behavior may be a dis- 
tinctive feature of minor mergers compared to secular processes 
in galaxy disks, which deserves to be modeled accurately in the 
near future. 

Minor mergers can quantitatively reproduce the increase in 
the rotational lag with height above and below the galaxy mid- 
plane that was fou nd for the Milky Way (IChiba & Beersll2000l: 
iGirard et ani2006h , but only if the satellite is accreted along a di- 
rect orbit; for retrograde orbits no trend of the rotational lag with 
z is found, up to distances of 5 kpc from the midplane, above 
which it even decreases slightly. Toget her with the increasin g 
trend of stellar eccentricities with height (iDi Matteo et al 
this is another piece of evidence suggesting that if the Milky Way 
thick disk has formed through heating of a pre-existing thin disk 
by minor merger(s), the orbit of the satellite(s) should have been 
prograde. 

Recently, iBournaud et al.l (l2009l) proposed a scenario in 
which stellar thick disks could have formed in gas-rich galaxies 
at high redshifts, through scattering of stars by massive clumps. 
Comparing the kinematics of minor merger remnants with those 
of thick disks formed in clumpy disks we find that, under certain 
conditions, both scenarios can reproduce a rotational lag of the 
old thin and old thick disk stars that are compatible with obser- 
vations for the Milky Way. In particular, disks heated by mas- 
sive clumps show a rotational lag which increases with height 
z, as found for minor mergers on direct orbits. Even so, the 
two processes leave distinctly different signatures in the verti- 
cal sur face profiles of the remnant disks, as shown in lOu et al.l 
(I2OIII) . Obviously, disks heated by massive clumps cannot lead 
to the counter-rotating thick d isks observed in external galaxies 
(lYoachim & Dalcantonll2008l) . For this to occur, an external ac- 
cretion of material (gas or stars) is necessary. 
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Fig. 19. Line-of- sight velocity, v, velocity dispersion, (Tv, and the 
v/(Tv ratio as a function of scaled radius in remnant stellar disks 
after both a single and two consecutive dissipative minor merg- 
ers. The same data are also shown for the initial stellar disks, 
before merging. The gas mass fraction of the merging galaxies 
increases from the left to the right, from for the SO case, to 0.2 
for the Sbs. 

Counter-rotating material is found in all our simulations 
where the satellite is accreted along a retrograde orbit. The 
amount of counter-rotating stars (all associated to stars initially 
in the satellite galaxy) increases with z as a result of the increas- 
ing contribution of satellite stars to the whole stellar content, and 
at z = 3 - 5 kpc the fraction of counter-rotating stars is between 
10% and 15%. We showed that the effect of multiple mergers 
on the generation of a counter-rotating disk is cumulative, in the 
sense that the fraction of counter-rotating stars doubles after the 
accretion of a second, comparatively massive, satellite on a retro- 
grade orbit. It is thus possible that some of the galaxies observed 
to ha ve a high fraction of counter -rotating material in the thick 
disk (lYoachim & Dalcantonl2005l) could have accreted a number 
of satellites on retrograde orbits. The fraction of counter-rotating 
stars we find in our minor merger simulations is quantitatively 
consistent with the amount of counter-rotation observed in the 
thick disks of nearby galaxies (lYoachim & Dalcantonll2005l) . 

Finally, we discussed the possibility of discriminating stars 
accreted from the satellite galaxy from those formed in the pri- 
m ary disk on the ba sis of their AM content, as recently proposed 
bv lVillalobos&Hel mi (2009). We find that in general the AM 
content is a better discriminant for stars located at vertical dis- 
tances from the galaxy midplane g reater than 1 kpc rather tha n 
at smaller heights, as proposed by IVillalobos & Helmil (l2009l) . 
Moreover, as the fraction of satellite stars increases with z in 
our models, their detection at high z should also be more likely. 
However, we want to point out that it may not always be possible 
to use the AM content to distinguish between stars originally in 
the primary and those accreted from the satellite because there 
are cases where no difference in their AM content is seen at any 
radius. 

Overall, the qualitative and some quantitative agreement 
with observations of the thick disk of the MW and other nearby 
galaxies suggest that minor merger remains a viable mechanism 
for forming the thick disk. With future observations, it may be 
possible to search for other signs of the impact of minor merg- 
ers on disk galaxies, such as determining the relative angular 
momentum of stars as a function of age in nearby disk galax- 
ies. Minor merger models, such as those discussed here, make 
specific predictions which can be tested observationally (e.g., by 
GAIA). 
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